circular convolution using matrix method